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Abstract 

A self-consistent particle-phonon coupling model is used to investigate the 
properties of the isobaric analog resonance in 208 Bi. It is shown that quanti- 
tative agreement with experimental data for the energy and the width can be 
obtained if the effects of isospin-breaking nuclear forces are included, in addi- 
tion to the Coulomb force effects. A connection between microscopic model 
predictions and doorway state approaches which make use of the isovector 
monopole resonance, is established via a phenomenological ansatz for the op- 
tical potential. 
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I. INTRODUCTION 



Since its discovery some 35 years ago, the isobaric analog resonance (IAR) has always 
attracted a considerable interest because it is one of the clearest manifestations of approx- 
imate isospin symmetry in atomic nuclei. The main causes of isospin symmetry breaking 
are essentially the long-range Coulomb force, and to a lesser extent the short-range nu- 
clear forces of the charge-symmetry breaking (CSB) type and charge-independence breaking 
(CIB) type PJ. Quoting a sentence from a review paper on this subject [[|, the Coulomb 
interaction "is strong enough so that its effects are quite easily detectable and weak enough 
not to destroy or alter considerably what the nuclear force has produced." Despite this 
favourable situation, we are still missing a microscopic model which is able to reproduce one 
of the physical observables related to the breaking of isospin symmetry in nuclei, namely 
the total width of the IAR, taking into account all the isospin-breaking effects. 

The general features of the IAR are qualitatively well understood. In the absence of 
isospin-breaking effects, the analog state would be degenerate with the ground state of the 
(N, Z) parent nucleus and it would have zero width. The effect of the Coulomb and other 
isospin-breaking forces is to shift up the analog state by several MeV (the so-called Coulomb 
displacement energy). However, these isospin-breaking forces do not induce strong isospin 
mixing of the states. The analog state, just as its parent state, has practically pure T = T 
isospin (To = (N — Z)/2) whereas the neighbouring states in the daughter nucleus have 
mostly T = Tq — 1. Thus, a simple argument based on the Fermi golden rule shows that a 
small spreading width for the IAR must be expected. In addition, the escape width must 
also be small because the main escape channels are those of low energy protons. 

When one attempts to build a microscopic model of IAR, the requirement that isospin 
symmetry must be restored if isospin-breaking forces are switched off should be taken into 
account. It is known that m N ^ Z systems the Hartree-Fock (HF) approximation in- 
troduces a spurious isospin symmetry breaking and that a self-consistent charge-exchange 
random phase approximation (RPA) restores this broken symmetry ||. By self-consistent 
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RPA we mean that the HF single-particle spectrum and the residual particle-hole interaction 
are derived from the same effective two-body force. Therefore, we shall consider here only 
this self-consistent framework. Charge-exchange RPA was applied for the first time to IAR 
studies in Ref. ||] (see also and references therein). Calculations of other types of exci- 
tations using charge-exchange RPA can be found in Ref. |J . However, the RPA description 
can at best give information on escape widths if continuum effects are included but it cannot 
shed any light on spreading widths because this spreading is due to states beyond the one 
particle-one hole (lp-lh) space. An extension of the model space to include 2p-2h configu- 
rations and leading to second RPA would be a more appropriate scheme. In heavy nuclei 
where the neutron excess is large one can safely replace the RPA by the Tamm-Dancoff 
approximation (TDA). In Ref. || such a second TDA calculation was performed for the 
IAR in 208 Bi and it gave a reasonable estimate of the escape width and spreading width 
T^. More recently, calculations based on a particle-phonon coupling model also led to 
satisfactory values of and but some discrepancies with experimental data still remain. 
In the above mentioned works, only isospin breaking due to the inclusion of the Coulomb 
force in the HF mean field was considered. 

The purpose of the present paper is twofold. Firstly, we wish to examine the effects 
of the isospin-breaking nuclear forces on the properties of the IAR. It is known that the 
CSB and CIB interactions affect the mass number dependence of the Coulomb displacement 
energy anomalies and bring them in better agreement with experiment || . Moreover, these 
isospin-breaking nuclear forces lead to an increase of about 50% of the isospin mixing in 
nuclear ground states ||. Thus, one should expect also some sizable effects on the values of 
Tj AR . We study these effects within the particle-phonon coupling model |7j together with 
the short range parametrizations of the CSB and CIB forces of Ref. ||. The calculations are 
performed for the typical case of 208 Bi. Secondly, it is instructive to study the connection 
between a microscopic model such as the present one or the second TDA model of Ref. || , 
and IVMR doorway state approaches [0,|n|] where the shift and width of the IAR result 
from the coupling via the Coulomb interaction of an ideal analog state with the isovector 
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monopole resonance (IVMR) playing the role of a doorway state. Here, we show that 
this connection can be established if one makes a phenomenological ansatz for the isospin 
dependence of the nucleon optical potential. 

In Sect. II the microscopic model is presented and its results are discussed in Sect. III. 



The connection between the microscopic model and the approach of Ref. [IT| is shown in 
Sect. IV. Conclusions are drawn in Sect. V. 

II. THE MICROSCOPIC MODEL 

The RPA extended so as to include the continuum coupling and the particle-phonon 
coupling has been described in detail in Ref. |7[] and therefore, we shall first recall only the 
main features of the model and then proceed to the specific points of the present calculation. 



Starting from an effective Hamiltonian H with two-body Skyrme interaction ||i2| , |T3 
Coulomb interaction and CSB-CIB interactions which will be described below, the HF 
equations determine the self-consistent mean field of the parent nucleus. This mean field 
is diagonalized on a basis of 15 shells of harmonic oscillator wave functions (hu osc = 6.2 
MeV for 208 Pb). This procedure provides a discrete set of levels (occupied and unoccupied). 
We select all occupied levels, and 6 unoccupied levels for each value of with increasing 
values of the radial quantum number n. This determines the subspace Q\ of discrete lp-lh 
(proton particle-neutron hole) configurations. 

To account for the escape width and spreading width of the IAR, we build two 
other subspaces P and Q2- The space P is made of particle- hole configurations where the 
particle is in an unbound state orthogonal to all the above discrete single-particle levels. The 
method to calculate these unbound states is described in Ref. 0. On the other hand, the 
space Q2 is built with the main configurations which are known to play a major role in the 
damping process of nuclear giant resonances: these configurations are lp-lh states coupled 
to a collective vibration. We have included in our model space all the isoscalar vibrations of 
multipolarity 2 + , 3~ and 4 + up to the energy of 20 MeV and which exhaust more than 1% 
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of the energy-weighted sum rule (EWSR) of the corresponding multipole operator. These 
collective vibrations are calculated consistently in HF-RPA within the Qi space and they 
constitute the phonons of our particle-phonon coupling model. 

Using the projection operator formalism one can easily find that the effects of coupling 
the subspaces P and Q 2 to Qi are described by the following effective Hamiltonian acting 
in the Q\ space: 

H (E) = QiHQx + W\E) + W\E) 

= Q 1 HQ 1 + Q,HP 1 PHQ, + Q,HQ 2 1 , Q 2 HQ 1 , 

E — PHP + ie E — Q2HQ2 + «e 

(1) 

where E is the excitation energy. For each value of E the RPA equations corresponding 
to this effective, complex Hamiltonian TC(E) are solved. The resulting sets of eigenstates 
enable one to calculate all relevant quantities such as giant resonances energies and widths 
(see Ref. [0]). In practice, we use the HF-RPA states (corresponding to positive and negative 
eigenvalues) as a basis for the Qi space because we can truncate this basis and neglect states 
which have negligible T_ strength. 

A simplifying approximation is made when calculating the matrix elements of by 
neglecting the interactions among the states within Q 2 . Each matrix element Wp hp , h , is a 
sum of four terms whose diagrammatic representation is shown in Fig. 1. To evaluate these 
diagrams we use the following expression for the particle-vibration vertices: 

V = Y. E (a\g^(r)v(r)Y LM (r)\/3) a\a p . (2) 

a/3 LnM 

This form comes from the particle-vibration coupling model where the vibration (phonon) 
\n) is characterized by its angular momentum L and its radial transition density Q^^r). 
The form factor v(r) appearing in (0) is related to the particle-hole interaction derived from 
the Skyrme force by Vphifx,?^) = v(r-\)8(fi — f^). When deriving this form factor, the 
velocity-dependent terms of the residual Skyrme interaction are dropped. 
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It was pointed out in Ref. that the particle-phonon coupling model does not auto- 
matically insure the isospin symmetry properties of the nuclear forces (contrarily to a fully 
microscopic 2p-2h model like in Ref. |J). One must therefore enforce isospin symmetry in 
the evaluation of Wp h ^ p , h , by an appropriate isospin projection procedure 0. 

In the present calculations the Skyrme interactions SIII |T2[ and SGII [[TJ| have been 



employed for the isospin symmetric part of the Hamiltonian. In addition, CSB and CIB 
effective nucleon-nucleon forces determined in Ref. || are also included. These forces were 
obtained from a short-range expansion of Yukawa-type interactions and they have a form 
similar to that of Skyrme forces: 

Vcsb = \ {r z (l) + r,(2)}{ So (l + y P a ) + + y x P„) 

x (k 2 + k' 2 ) + s 2 (l + y 2 P cr )k'-k}, (3) 

and 

Vcib = \ t z {1)t z {2){u {1 + z Pa) + \u x {\ + Zl P a ) 

x(k 2 + k' 2 )+u 2 (l + z 2 P (T )k'-k}. (4) 

The parameters s$ and are given in Ref. P|, and all exchange parameters j/j and Zi are -1 
because of the singlet-even character of Vcsb and Vcib- Therefore, they do not contribute as 
residual particle-hole interactions in the isovector channel and their only influence is through 
their contributions to the HF mean field. 



III. DISCUSSION OF RESULTS 

Three types of calculations have been done which are labeled by I, II, III respectively. 
In calculation I the starting Hamiltonian contains only the Skyrme interaction without the 
Coulomb force between protons and without CSB-CIB nuclear forces. In calculation II the 
Coulomb interaction is also included. In calculation III the CSB-CIB forces are added. Thus, 
the three calculations have an increasing degree of isospin breaking and they are expected 
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to lead to increasing values of the IAR width according to the arguments recalled before. 
In calculations II and III the Coulomb exchange contributions to the mean field are treated 
in the Slater approximation whereas the Coulomb p-h residual interaction is dropped. The 
results obtained by using respectively the interactions SIII and SGII are shown in Tables 1 
and 2. Calculation II is in principle equivalent to what was done in Ref. J7| except for two 
changes in the model spaces P and Q2, and in the averaging parameter A. Because of these 
changes, the values of the mean energies and widths of the IAR we quote as calculation II 
are slightly different from the correspnding values reported in . The spaces P and Q2 have 
been enlarged with respect to the calculation of Ref. by including a larger (96 instead 
of 64) number of unoccupied proton single-particle states. This affects mainly the escape 
and spreading widths calculated with the interaction SIII which both increase by about 10 
keV whereas the widths calculated with the interaction SGII are less affected. Also, the 
averaging parameter A (see Ref. j7|) employed in the calculation has been changed in the 
present work from 100 keV, as it was in 0, to 200 keV. In fact, we have carefully studied the 
dependence of the spreading width on the parameter A, and we have found that IN- increases 
by about 20 keV when A is changed from 100 keV to 200 keV (as it was already discussed 
in Ref. |7|) and also increases by another 20 keV when A is changed from 200 keV to 400 
keV. On the other hand, 1^ remains constant if A is set above 400 keV. Therefore, we have 
adopted the value of 200 keV for A since this value is intermediate between keV and 400 
keV (the value at which saturates). This study allows us to say that the uncertainity on 
the values of 1^ due to the freedom of choice of the parameter A is about ± 20 keV. 

The calculation without any isospin-breaking force (calculation I) has been performed to 
show that our procedure can recover the isospin symmetry reasonably well at each step. The 
first step (a) is discrete TDA since discrete RPA is not possible in this case because of the 
negative energy configurations. Indeed, without the Coulomb interaction the neutron states 
lie higher in energy than their proton partners and the energies e p — Eh of the unperturbed 
proton particle-neutron hole + configurations are negative. For the six excess neutron 
hole levels and the corresponding proton particle levels, these values are between -7.3 MeV 
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and -8.3 MeV in the case of the force SGII and between -8.6 MeV and -10 MeV in the 
case of the force SIII. In an ideally accurate numerical calculation these negative energy 
configurations would be coherently pushed up to form a collective state at zero energy. 
Small numerical inaccuracies can weaken the residual interaction, and as a result complex 
RPA eigenvalues may appear. This can be easily understood in the simple case of a schematic 
model, or in a case containing a single state, where the RPA and TDA energies are related 
by E 2 RPA = Ej> DA — V 2 with V representing the interaction term. In the limit Vc = 
the r.h.s. may become slightly negative and consequently Erpa is imaginary, if E\ DA and 
V 2 do not exactly cancel numerically. This difficulty of complex RPA solutions in the case 
Vc = is well known and, for example, in Ref. || the way out was to slightly renormalize 
the nuclear residual interaction. Therefore, only TDA is possible for the study of the IAR 
without the Coulomb interaction. One should also note that, contrary to a naive view, 
adding the Coulomb interaction does not result in an overall shift of the RPA eigenvalues. 
This can again be seen, e.g., in the framework of a schematic model. The second step (b) is 
RPA with the coupling to the continuum, which gives no width since the IAR is below the 
proton emission threshold. The aim of the third step (c) (inclusion of the spreading width) 
is to know whether the particle-phonon coupling model we have adopted can introduce 
some spurious width because of the approximations made: finite size of the set of lp-lh 
plus phonon states, overcompleteness of the lp-lh plus phonon basis, violation of the Pauli 
principle, simple form of the particle-phonon vertex function. The fact that we obtain only 
4 keV for the spurious ]> with the force SGII is very satisfactory since it means that the 
approximations we have mentioned are safe in this case. 

In the case of the force SIII the spurious width (24 keV) is larger than that obtained 
with the force SGII. We have reached the conclusion that this spuriosity is due by half 
to the particle-phonon coupling and by half to the fact that HF-TDA is not completely 
able to recover the isospin symmetry. To estimate the spuriosity due to the particle-phonon 
coupling, we have considered an ideal analog state at zero energy with the following schematic 
wave function: 
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1^4) 



1 



f_\HF)= X ; 




(5) 



with v 1 restricted to the excess neutrons and 



X 



1 



S(l v ,l v )S(J v ,j v ) y/2j„ + 1 / dr u n (r)u u (r) 



(6) 



We have calculated the coupling of this schematic IAR with the lp-lh plus phonon states 
adopted for all the calculations of Table I, and we have obtained a state at 27.6 keV whose 
width is 13 keV. Therefore, 13 keV is the broadening of the IAR introduced spuriously by the 
coupling with particle- hole-phono n configurations. The remaining 11 keV is still larger than 
the value obtained with the interaction SGII. This is on one hand related to the different 
single-particle levels obtained for the two interactions. As mentioned above, proton particle 
and neutron hole levels are more separated in the case of the force SIII. Consequently, the 
isospin breaking in the HF field is restored less efficiently by TDA and the energy of the 
IAR is less close to zero than in the case of the force SGII. Therefore, in the case of SIII we 
diagonalize the effective Hamiltonian (U) at a higher value of the energy E and the spreading 
width of the IAR at higher excitation energy is larger. Moreover, the imaginary parts of 
the self-energy terms in our model would cancel exactly to give a zero spreading width if 
the single-particle radial wave functions of neutrons and protons with the same quantum 
numbers (n, l,j) were identical. This cancellation can be seen by looking at the expressions 
of the four diagrams of Fig. 1 . However, the cancellation is not complete if the single- 
particle radial wave functions are different. The difference between radial wave functions is 
larger in the case of the force SIII, so again the spurious width is expected to be larger than 
that of SGII. 

If we now include the Coulomb interaction between protons in the HF mean field, the 
proton levels are pushed up and become higher than the corresponding neutron levels. The 
energies of the unperturbed proton particle-neutron hole + configurations have positive 
values between 11 and 11.8 MeV for SGII and between 8.8 and 10.3 MeV for SIII (for the 
six main configurations already considered above). This difference between the two forces is 
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essentially related to the fact that in the case of the force SGII the neutron holes are more 
bound. Therefore, the proton particles which enter the IAR wave function have less energy 
available (e p = Ejar + Sh) and a smaller probability to escape. This explains why the is 
considerably smaller than in the case of the force SIII. 

Finally, the most important result of our calculation is that the total width obtained by 
employing the force SIII (in the case of the complete calculation III, last column of Table 1) 
nicely agrees with the experimental finding T^or = 230 keV. The improvement with respect 
to calculation II (without CSB-CIB forces) is about 15%. This shows that CSB-CIB forces 
can contribute significantly to the total width of the IAR. 



IV. COMPARISON WITH THE IVMR DOORWAY STATE APPROACH 

We have seen that microscopic approaches, like the particle-phonon coupling model de- 
scribed above or the 2p-2h TDA model of Ref. || can give a reasonable description of IAR 
widths in spite of some sensitivity to the effective interactions. The spreading widths come 
from the coupling terms Q\HQ 2 appearing in of Eq. ([!]), and this coupling between the 
simple Qi configurations and more complex Q 2 configurations is produced mostly by the 
isosp in-conserving Skyrme interaction (in fact, this is the only residual interaction we keep 
in the calculation of in the previous sections). Thus, the fact that the resulting 1^ is 
non-zero is due entirely to the effects of the Coulomb and other isospin-breaking forces in the 
mean field as they produce a finite density of states with the isospin of the parent nucleus at 
the IAR energy. On the other hand, in the approaches of Refs. []T0|,p|,pT| originates from 



the coupling of an ideal analog state \A) with a specific doorway state, namely the IVMR in 
the daughter nucleus, via the isospin-breaking part of the Hamiltonian (usually the isovector 
component of the Coulomb force). Here, we show that these apparently different points of 
view can be connected. 

Let us first recall the expression for 1^ obtained in Ref. |JL1||. The Hamiltonian is assumed 
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to be a sum of an isospin-conserving part plus the Coulomb interaction, H = Hq +Vo, and 
the parent ground state which is eigenstate of H is denoted by |0). The three isospin 
components of the IVMR in the daughter nucleus are schematically written as 

\M;T + i,T -l) = |{|0) T=T »® Iph-Y^}^} , (7) 

where i = —1, 0, 1 and \ph~ l ) stands for a combination of monopole p-h excitations. In Ref. 
|TT| the IAR spreading width (we use this notation for the spreading width calculated 
by following the doorway state approach) was expressed in terms of the analog state energy 
E A , IVMR energies E^} + \ the width T M (E A ) of the IVMR evaluated at energy E = E A , 
and the reduced Coulomb matrix element: 

vc = ^ h ~ l ) T=X II V c ] II 0) , (8) 

where Vq is the isovector part of the Coulomb potential. If one neglects the isospin splittings 
of the IVMR and adopts a common value E^ +l ~ Em the expression of f ^ takes the simple 
form 

Furthermore, it was shown that the isospin mixing probability of the Tq + 1 component of 
IVMR in the parent ground state \ir) is given in second order perturbation theory by 

|cTo+l|2 "2(^Ti) (7r|T - T+|7r) 

(10) 

where AEm is the excitation energy of the IVMR in the parent nucleus which is approxi- 
mately equal to E A — E M of Eq. (|D|). We can safely neglect T M in Eq. (|S|) and thus obtain 



1 


vc\ 


2 


T + l 


AE M \ 


2 



, - lr M (E A )(n\T_T + \n) . (11) 



'4 



2 

In the microscopic models |||7| the spreading width of the IAR results from the couplings 
mediated by the isospin-conserving operator defined in section II. Denoting the RPA 
eigenstate corresponding to the IAR by \A), we can write its width as 
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= -2 Im (A\W l \A). (12) 



The IAR wave function can be well approximated by 



For the isospin-conserving interaction we make the ansatz 



\A) = -U= r-K) . (is) 



W l = a(f -f -b) 

= (a R + i ai )(f-f-b) , (14) 



and obtain 



(A\W'\A) = ^(n\T + (f.T-b)T_\n) 



2T f 



o 



±^\T + T-{f ■ f - b)\n) 

^(7r\T + T_\n}(n\(f.f-b)\n) . (15) 



o 



In the last step of the above equation the contributions of the excited states to the closure 
relation have been dropped because they are of order |ct +i| 2 - Let us introduce T by: 

(ir\f ■ f |tt) = f (f + 1) , (16) 

(T differs slightly from T because \ir) has isospin mixing) and choose b = T 2 . Then: 

(A\W l \A) = a[f{f + l)-b}=af. (17) 

Thus, the two expressions ( pT| ) and ([12|) are equal if the following condition is satisfied: 

a/ = -^^<7r|T_T + K>. (18) 

Next, we diagonalize the TDA schematic model with the interaction 

v eff = {a R + ia I ){f-f-b). (19) 

The complex eigenvalue is 

12 



Ea - i V -f = e ph + V -f- EK(P^" 1 ) T=1 ' T * = " 1 l r -l°)l 2 

ph 

= £ph + ^--2T + i ^.2T , (20) 

where e p h is the degenerate unperturbed energy of the + proton particle-neutron hole 
configurations. Replacing a/ by its value flT8|) we see that Ta becomes 



T A = h M (E A )^(n\f_f + \n) , (21) 

which is consistent with Eq. (jTT|) . Thus, the adopted interaction results in a value of 
the spreading width coming from the coupling of the IAR with states of 2p-2h type through 
the nuclear interaction, which is comparable with the result of the IVMR doorway state 
approach in which the spreading width is obtained through the coupling of the IAR to the 
IVMR due to the Coulomb interaction. It is left as a future problem the justification of the 
choice of (O) from a microscopic study of optical potentials. 



V. CONCLUSION 

Within the framework of a microscopic model based on self-consistent HF-RPA plus 
coupling with continuum configurations as well as with lp-lh plus phonon configurations, 
we have calculated the total width of the IAR in 208 Bi. We have shown that if the nuclear 
isospin-breaking forces of CSB and CIB type are included in the Hamiltonian in addition to 
the Coulomb interaction, the width of the IAR is increased by 15-20%. Thus, the nuclear 
isospin-breaking interactions which were already known to increase the isospin mixing of 
ground states have also significant contributions to the total width of the IAR. As far as 
comparison with experiment is concerned, the values of 1^ and Ytot calculated with SIII 
are in satisfactory agreement with the data whereas those obtained with SGII are not so 
good due to the peculiarities of the single-particle spectra of SGII. 

Our microscopic model introduces some spuriosity in the evaluation of the total width of 
the IAR. This spuriosity turns out to be quite small and we have shown that it is due partly to 
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the incomplete restoration of symmetry by TDA and partly to the lp-lh plus phonon model. 
Indeed, this model does not have the full self-consistency of a second RPA calculation which 
would be free in principle of spurious isospin violations. However, a second RPA calculation 
including also continuum effects would be extremely difficult and it has never been done so 
far. 

Finally, we have been able to propose for the first time a connection between the micro- 
scopic model and IVMR doorway state approaches for the spreading width of the IAR. This 
connection is possible by making a phenomenological ansatz for the isospin dependence of 
the nucleon optical potential. 
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FIGURES 

FIG. 1. Diagrams corresponding to the coupling of the lp-lh configurations to the more 
complicated states including a phonon (wavy line). The sum of the four diagrams gives the matrix 
element , h , quoted in the text. 
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TABLES 

TABLE I. IAR results with the interaction SIII. Three different types of interactions based on 
SIII are used in the calculations: I) Skyrme interaction without Coulomb force; II) Skyrme interac- 
tion with Coulomb force; III) Skyrme interaction with Coulomb force and CSB-CIB forces. Three 
different microscopic models are also adopted: a) TDA without the coupling to the continuum; b) 
RPA with the coupling to the continuum; c) RPA with the couplings to both the continuum and 
the phonons. Energies are given in MeV and widths are in keV. The percentage of total strength 
nio=(N-Z)/2 exhausted by the IAR is also shown. 
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TABLE II. IAR results obtained with the interaction SGII. For details, see the caption to the 
previous table. 



(a) Discrete TDA (b) RPA + W T (c) RPA + W T + W l 
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